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We introduce the minimal momentum subtraction (MiniMOM) scheme for QCD. Its definition allows the strong coupling to be 
fixed solely through a determination of the gluon and ghost propagators. In Landau gauge this scheme has been implicit in the early 
studies of these propagators, especially in relation to their non-perturbative behaviour in the infrared and the associated infrared 
fixed-point. Here we concentrate on its perturbative use. We give the explicit perturbative definition of the scheme and the relation 
of its /^-function and running coupling to the MS scheme up to 4-loop order in general covariant gauges. We also demonstrate, by 
considering a selection of Nj = 3 examples, that the apparent convergence of the relevant perturbative series can in some (though 
not all) cases be significantly improved by re-expanding the MS coupling version of this series in terms of the MiniMOM coupling, 
making the MiniMOM coupling also of potential interest in certain phenomenological applications. 

Key words: strong coupling constant, 4-loop running, minimal subtraction, momentum subtraction 
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1. Introduction 

The coupling constant of the strong interaction, a s , is one of 
the fundamental parameters of the Standard Model. It is spec- 
ified by its value in a particular renormalisation scheme at a 
chosen reference scale, p, conventionally taken to be the MS 
scheme at the Nf = 5 reference scale p - Mz- A recent assess- 
ment of experimental results yields af S (M Z ) = 0.1189(10) lUl]. 
This result, which is little changed if more recent experimen- 
tal results are taken into account (see, e.g., |2[] and references 
therein), is in excellent agreement with two recent, slightly dif- 
ferent, lattice determinations [2, 3] based on lattice perturbation 
theory analyses of short-distance-sensitive lattice observables 
computed using the MILC Nf = 2 + 1 configurations. 

Other schemes than the MS scheme are of course also pos- 
sible. For example it has been proposed in [4] that a particular 
product of dimensionless gluon and ghost dressing functions, 
Z and G, in the Landau gauge can be used to define a non- 
perturbative running coupling via 



«rv) = f-z(/)GV) 



(1) 



where g 2 = g 2 (p) is the renormalised coupling at scale p and 
the renormalisation condition, 



Z{p 2 )G\p% =ii2 = 1 , 



(2) 



is assumed. The definition (0} is particularly convenient since 
it allows the coupling and hence A MS to be determined by mea- 
suring two-point functions on a lattice. Our first steps towards 
such a determination for Nf — 0, 2 were presented at the 2007 



Lattice Conference |5J|. Further updates of these preliminary 
results were reported last year JaLZI]. This project, which is on- 
going, has the potential to provide an independent precision de- 
termination of a s from lattice simulations at purely perturbative 
scales, typically between 10 and 100 GeV JHQ As an important 
supplement we specify here the details of the renormalisation 
scheme underlying the coupling (HJ for QCD in general covari- 
ant gauges. We also provide the explicit 4-loop /3 function for 
o/f M in all such gauges. This information will be important to 
the previously mentioned lattice analysis and to quantifying the 
truncation error on the resulting a s determination. 

The product in (HJ is dimensionless and renormalisation 
group invariant, and it reduces to the running coupling of a per- 
turbative momentum subtraction scheme (MOM) at large p 2 . 
This makes it a suitable candidate for a non-perturbative exten- 
sion, though such extensions are, of course, not unique. The 
underlying renormalisation condition (f2]i respects infrared scal- 
ing (with 0.5 < k < 1 H), 



Zip 2 ) 
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and G(p 2 ) 



(p 2 /Kc D r 



(3) 



for p 2 — » 0, as predicted by a variety of functional con- 
tinuum methods for Landau gauge QCD including studies of 
Dyson-Schwinger Equations (DSEs) Stochastic Quantisa- 
tion [ 10], and of the Functional Renormalisation Group Equa- 
tions (FRGEs) ifTlll . This conformal infrared behaviour of the 



'The QCD scale parameter of the underlying scheme (the MiniMOM 
scheme, see below) is roughly 450 MeV for Nf = 0, or 430 MeV for Nf = 2. 
Non-perturbative contributions to the gluon and ghost dressing functions are at 
least of the order (Nqco/p) 2 - They are negligible at scales ft above 10 GeV. 
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purely gluonic correlations in Landau gauge QCD is consistent 
with the conditions for confinement in local quantum field the- 
ory I HI 131. but it is yet to be observed in lattice simulations^ 

If the infrared scaling behaviour ([3]) is realised, the renormal- 
isation condition (|2]i holds beyond perturbation theory and can 
be imposed at any (space-like) subtraction point p 2 — p 2 > 0. 
This is because the running coupling defined by flj then ap- 
proaches a finite infrared fixed-point, a™ M — » a c for p 1 — > 0, 
with a c a 8.9 /N c obtained under a mild regularity assumption 
on the ghost-gluon vertex ^ . For the purposes of this Letter it 
suffices, however, that the coupling (H}, is well-defined pertur- 
batively, independent of the infrared scaling behaviour in (01. 

A special feature of Landau gauge which underlies the def- 
inition (Q3 is the non-renormalisation of the ghost-gluon ver- 
tex ||l5ll . The possibility of taking advantage of this non- 
renormalisation has been criticised in the past on the 
grounds that the Landau-gauge ghost-gluon vertex acquires a 
momentum dependence in all common MOM schemes, already 
at one-loop level 11811 . despite the absence of ultraviolet diver- 
gences in this vertex in Landau gauge. There is, however, no 
contradiction here at all, as will be explained in the next section. 
The basic reason is that the non-renormalisation of the ghost- 
gluon vertex ensures the existence of a scheme for which the 
notion of ([TJ as a running coupling makes sense and which this 
running coupling implicitly defines, without the need to use an 
asymmetric momentum scheme fl^l . We call it the MiniMOM 
scheme for reasons that will become clear below. A useful fea- 
ture of the MiniMOM scheme is that it can be related to the 
MS scheme at four-loops Ja, |7[] without the need to compute 
vertices to three loops in perturbation theory. 

After providing more details on the MiniMOM coupling and 
its relation to the MS coupling below, we consider, in Sec. [4] 
the MS and MiniMOM versions of the perturbative series enter- 
ing a selection of phenomenological applications, demonstrat- 
ing that, in some cases, the apparent convergence of the series 
is much improved by the use of the MiniMOM coupling. Our 
conclusions and outlook are given in Sec. [5] 

2. The minimal MOM scheme 

Some confusion in the literature concerning the running cou- 
pling of Eq. ([D arose in relation to the misconception that this 
definition rests on the non-renormalisation of the ghost-gluon 
vertex in Landau gauge |15|]. This definition is not, however, 
based on requiring that the ghost-gluon vertex reduce to the 
tree-level one at a symmetric subtraction point k 2 — p 2 = q 2 = 
p 2 . In particular, the scheme underlying (Q} is not the MOMh 
scheme of Ref. [20], but is, instead, a minimal MOM (Mini- 
MOM) scheme, which is defined as follows: 

As in every MOM scheme, the gluon and ghost renormalisa- 
tion constants Z3 and Z3 are defined by requiring 



the perturbative realisation of (f2j) valid for p? » A 2 where A 
is the scale parameter of the scheme. Z(p 2 ) and G(p 2 ) are the 
dressing functions in the renormalised gluon and ghost prop- 
agators of Landau-gauge QCD, which are in (Euclidean) mo- 
mentum space of the form 



p li p v \ Z(p 2 ) 



and 



D' 



ab 



a(P) = S 



ab 



G(p 2 ) 



(5) 



(6) 



respectively. Instead of imposing additional analogous renor- 
malisation conditions on vertex functions, requiring certain ver- 
tex structures to equal their tree-level counter parts at some 
symmetric or asymmetric subtraction point, we here supple- 
ment (|4|l by the further condition 



^1 



(7) 



where Z\ is the renormalisation constant of the ghost-gluon ver- 
tex, whose momentum dependence is thus the same as in the 
minimal subtraction schemes. The renormalisation constants 
for the three and four gluon vertex functions are then deter- 
mined by the Slavnov-Taylor identities as usual, 



Z3 — Z3 — 5 

Zi = —Z\ and Z4 = — Zf 

z, z? 



(8) 



The extension to quarks is straightforward, and this defines 
the MiniMOM scheme: momentum subtraction for the gluon, 
ghost and quark propagators, and minimal (MS) subtraction 
for the ghost-gluon vertex together with the Slavnov-Taylor 
identities to fix the remaining vertices, including the quark- 
gluon vertex. Note that the renormalisation constants for these 
remaining vertices then differ from those in the MS scheme by 
ratios of the propagator (field) renormalisation constants in the 
MiniMOM (MM) and MS schemes, for example, 



7 MM 7MS 

^3 _2 Z MS 

^MS 2 MM ' 



7MM 

and z 4 MM = ^= 

^3 



7MS \ 



7' 



MM 



7-MS 



(9) 

Likewise, from the general definition of the renormalised cou- 
pling constant, 



a s (p) = 



An 



Z3Z 2 
Z 2 



'bare 

An 



(10) 



the MiniMOM and MS scheme couplings, with the definition 
Z™ = zfs in O, are related by 



Z(p 2 



1 and G(p 2 



1 



(4) 



2 In order to do that one needs a proper non-perturbative definition of BRST 
symmetry on the lattice which is possible in principle, but not realised in present 
lattice implementations of Landau gauge. For a recent discussion, see I14ll . 
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Note that for this conversion we only need to know the pertur- 
bative expansions of the ghost and gluon propagators but not of 



2 



any vertex structures. Rather, with the MOM renormalisation 
conditions for these propagators or their non-perturbative 
extension (f2]l for that matter, we simply obtain, 



V) 



„MS 



(M) 



Z(^ 2 ) MS G 2 (^ 2 ) MS 



(12) 



from the gluon and ghost dressing functions evaluated at p 2 = 
jj 2 in the MS scheme. Alternatively, we can relate the Mini- 
MOM coupling to the MOMh scheme just as easily via 



(Y 



MM 



0") = 



^MOMh ^ 
"1 



7 MS 



af 0Mh (fi) . 



(13) 



All these conversion identities are valid for arbitrary linear- 
covariant gauges and not restricted to Landau gauge. The spe- 
cial feature of Landau gauge is that there Z^ s = Z^ M = 1. Be- 
cause the Landau gauge ghost-gluon vertex trivially reduces to 
its tree-level form when one of the ghost momenta is set to zero, 
the MiniMOM scheme then agrees with the asymmetric MOMh 
scheme of Ref. [21] (called the T scheme in Ref. [19]) which is 
defined from renormalising precisely this vertex structure. The 
MiniMOM scheme is defined, however, so as to not require 
knowledge of any vertex structure beyond the MS scheme con- 
tributions as determined by their ultra-violet divergences. This 
makes it particularly useful for a lattice determination of a s 
from the perturbative behaviour of QCD Green's functions as 
described in Refs. JUSH. 

For early 2 and 3-loop calculations in Feynman gauge, see 
1 23]. The complete 2-loop results for general gauge parameters 
are given in 12411 . To obtain the 4-loop version of the conversion 
between a^ s and o'^ IM from Eq. (fT2l . we use the 3-loop expres- 
sions for the gluon and ghost self-energies found in Appendix C 
of 121'] which, with a = af s /n, for Nf massless quarks yields 

a MM /Q/ MS _ j + Di a + Di a 2 + ^ a 3 + 0^4) > (U) 

where, with ( 3 = 1.2020569, & = 1.0369278, and C F = (N 2 - 

l)/2N c , 



= flio + d n N f , 
. = d 2 o + d 2 \Nf + d 22 N 2 , 

j _ [76063 39 r ,( 269 
"20 - [50736 - ' 1 

(M 



d u - 



dio = 
D 2 = d 2Q 

j _ [ 76063 39 r , I 269 , , 
"20 - [50736 _ 128 « + (, 2304 + 646) £ + 

(256 _ T28^ 3 )^ + 256 ^]^c> 

d 2 \ = -[I - k&] c f ~ [M + \& + m £\ N c < 
di2 = H . 
D 3 = J 30 + dnNf + di 2 Nj + d^Nj , 

j _ [ 42074947 _ 20225 r _ 7805 y , 
"30 [ 2985984 9216^3 12288^ 5 " 1 " 

I 17743 ' 



(15a) 



(15b) 



(15c) 



17743 , 10335 r _ 221/ W-l 

11472 + 9216 ^3 1024^> 5 /S + 

/ 16235 _ _71_ r _ YVS_ r \ A , 

\ 36864 1536^3 6144^? + 



Table 1 : The numerical values of the in the 4-loop relation between af M 
and Landau gauge, Eqs. 1141 and 115\ with f = 0, for N c = 3 and 

Nf = 0, 2, and 3. 





N f = 


Nf = 2 


N f = 3 




3.52083333 


2.96527778 


2.6875 


D 2 


29.71718945 


20.96900712 


16.82639745 


D 3 


291.4436449 


175.9308786 


127.6686773 



/ 1207 
^ 12288 

/ 169 
^ 12288 



3072^3 3072 ^ 5 )^ + 

\ 12288 ~~ 3072^3 + 72288 ^J^c 

j _ [ _ 202997 _ 217 r , 5_r _ 
"31 _ [ 31104 288^3 f \ 2 <,5 



+ »f3)f-(5^-A6)^K- 
...... Bft -f?5 + (m-T 6 b)f]CfN c < 

[ 143 , 37 r 5 r l r 2 
[576 + 48« - j(5\C f , 

*2 = [Hi + fa£s + (5^2 + ]w c + 

r 7nni n r ] r j _ 125 

[«JW . "33 - -5832 ■ 



I 505 
^5184 

[ 41243 
[ 10368 

[ 143 



[ 7001 _ 23 j 
[ 10368 24- 



This conversion depends on the gauge parameter £ = ^ as in 
every other momentum subtraction scheme, though in the Mini- 
MOM scheme, this dependence is comparatively weak. At one- 
loop level, for example, the coefficient of the leading £ depen- 
dence around £ = in the above conversion is 3 times smaller 
than that of the asymmetric MOMh scheme (which coincides 
with the MiniMOM scheme at £ = 0). 

For f = the same conversion can be obtained from the 
product of the scheme-invariant propagators given for Landau 
gauge in Sec. 4 of [25]. We verified the general result $15[ for 
N c — 3 and Nf = 0, 3 and 6 from the explicit expressions given 
there. Another check of our conversion, up to and including 
0(a 2 ), can be obtained using Eq. ( fT3l in Landau gauge, where 
= 1, with the 2-loop expression for_Z^ OMh (fi) = F ^ 2 ) 
in the symmetric MOMh scheme from 12011 . together with the 
3-loop conversion from a^" s to a^ I0Mh as given there 

The numerical values of the Dk in Landau gauge are given 
explicitly, for N c = 3 and Nf = 0, 2, and 3, in Tab. 03 As an 
illustration of the weakness of the dependence on the corre- 
sponding N c = 3, Nf = 3 results are 

D\ = 2.6875 + 0.3750£ + 0.1 875£ 2 , (16) 

D 2 = 16.8264 + 2.5977^ + 0.9769£ 2 + 0.1758£ 3 , 

D 3 = 127.6687 + 26.7739^+ 8.3907£ 2 + 1.9621^ 3 +0.3349^, 

with similar, slightly less ^-dependent, results for Nf = 0, 2. 

In order to compare actual values of a s in the MiniMOM 
scheme (at £ = 0) to the corresponding ones in the MS scheme 
we give two important examples. First, at the mass of the Z bo- 
son, m z = 91.2 GeV with op(m|) = 0.1189 H, from Eq. 



3 There appears to be a typo in Table 2 of 1 20], the entry for the 2-loop 
coefficient of T proportional to N c Nf should probably read -0.115(2) instead 
of -0.151(2). 



we obtain 



af M (m|) = 1.096af s (ml), 



(17) 



for Nf = 5, while at the mass of the t lepton, m T = 1.777 GeV, 



af M ('«?)= l.59af s (m 2 T ). 



(18) 



where, to be specific, we have used of (m 2 ) — 0.322, the value 



obtained by running (m z z ) = 0.1189 down to the Nf = 3 
scale m 2 using the standard 4-loop running fl26ll . 

We conclude this section with a few comments on quark 
mass effects. As in every off-shell subtraction scheme, the run- 
ning coupling and beta function of the MiniMOM scheme in 
principle depend on the quark masses. This is evident from 
its relation to the (mass independent) MS scheme, Eq. ( fT2l . 
in which the MS scheme gluon and ghost dressing functions 
will depend on the masses in the quark loops. These have not 
been included and our conversion formulas are therefore strictly 
speaking valid only for Nf massless quarks. To fully account 
for finite quark masses at this level one would need the corre- 
sponding 3-loop expressions for the gluon and ghost dressing 
functions with massive quark loops which have not been calcu- 
lated to our knowledge as yet. 

We can estimate the leading quark mass effects, however, 
which will affect the conversion formula at the 2-loop level. 
These are obtained from the 1-loop vacuum polarisation with 
massive fermions in the gluon self-energy, and they lead to an 
increase of D\ (via dn) as compared to the massless case. For 
each quark flavour with mass m q one then separately obtains a 
transcendental function d u (y) of y = m 2 /^ 2 which approaches 
-5/18 fory — > 0, i.e., for fi 2 » m 2 . Using the current upper 
limits from the Particle Data Group for the average up/down 
mass of 5 MeV and the strange mass of 130 MeV as commonly 
given at // = 2 GeV, we observe a maximum increase in D\ by 
0.15% at n — 2 GeV as compared to the massless Nf — 3 value 
given in Tab. Q] At fx = 1 GeV, with correspondingly larger 
light quark masses, the same upper bound for the increase in 
D[ reaches 1%. At larger scales the effect rapidly decreases. 
In particular, the explicit comparisons in Eqs. dT7T > and ( fT8l re- 
main unaffected by the corresponding changes in D\ (note that 
even with the charm and bottom quark masses included, the in- 
crease in D\ will be less than 0.1% at ix = mz as compared to 
the Nf — 5 massless flavour value there). 

The most noticeable quark mass effects will of course occur 
right at the decoupling scales /j = m q . At the charm threshold, 
for example, with jx = m c = 1 .27 GeV, the quark mass contri- 
butions to the vacuum polarisation lead to an increase in D\ by 
around 13% as compared to the massless Nf - 4 value (at fi — 2 
GeV this increase is reduced by a factor of 2 already). Consid- 
erable charm-quark mass effects should therefore be expected 
when converting the 4 flavour MiniMOM coupling to the MS 
coupling in the phenomenoligically interesting range between 
1 and 4 GeV. Without a more detailed knowledge of these ef- 
fects an Nf -4 conversion would therefore not be advisable. 
Here we use the MiniMOM to MS conversion up to at most 
Nf = 3 for which the quark mass effects are very small. For 
now, matching to the Nf — 4 and 5 regimes should always be 
done after the conversion, for the MS coupling in the usual way. 



There are no charm and bottom quarks in the lattice deter- 
minations of the MiniMOM coupling which are presently re- 
stricted to Nf — and 2 and which will be extended to 2 + 1 
light flavours in due course. At the relevant high momentum 
scales quark mass effects should then be completely negligible 
in the conversion to a^. In addition, it is always a possibility 
to remove any small residual light-quark mass effects by extrap- 
olation, if necessary. 

3. Beta-function coefficients of the MiniMOM coupling 

The running of the coupling constant as the scale, [i, changes 
is controlled by the (scheme-dependent) B function which, at 
small couplings, is defined by 

2 da(fi 2 ) 



d/u 2 



= t3{a) := 



i+2 



(19) 



where a = a s /n. The B function in the MS scheme is known 
to 4-loop order, the expressions for the corresponding coeffi- 
cients B^ s , k — 0, • • • , 3, for general N c and Nf, being given 



in Refs. 12711 . These results, together with the relation between 
the MS and MiniMOM coupling given in Eqs. dT4b and ( IT3T >. 
and the 3-loop version of the expression /u 2 d^/dfi 2 - ijyj, for 
the running of the renormalised MS gauge parameter, £ with 
73 the gluon anomalous dimension, the 3-loop expression for 
which can be found in Appendix D of Ref. 12111 . allow us to 
obtain the B function coefficients in the MiniMOM scheme to 
4-loop order. For general N c and Nf, we find 



fl MM 

Po 



B^-lNf] 



(20a) 



< M = I [t N c ~ lNfN c - C f Nf] + Bf + B^Nf , (20b) 

oMM _ [9655 143 y 1 w 3 [2009 , 137 y 1 A j AT 2 

Pi ~ [4608 - 512 ^ 3 J N c ~ [2304 + 768 «J NfN c (20c) 

+ [m + T4 ft] N 2 f N c £,] CfN f N c 

+ a^/ + [I-56]^ 



oMM _ [ 1381429 _ 225335 s- _ 85855 p ] at4 
^3 [ 165888 110592 ^ 3 73728 ^5 J 

[ 244549 , 3395 r , 35965 /■ ] M N 3 
^ [ 55296 18432 ^ 3 36864 "J ly f ly c 

• [ 5 , 11 y 1 60685 85 y 55 y \ r A7 
+ [ _ 96 + T « _ [ 18432 - 64 ft _ 48 ft j C f N f 

B-ilftK+(j§ + f ft-ift)^ 

, / 2357 43 y 5 y \ r at2 
+ 12304 ~ 96 ft _ 2A^) L f N f 

-(m + mb) N f] N e 



(20d) 



V648 ' 432 * 3 J"f\ ' c 

+ W 6 C f Nf + {^-±b)N 2 f 

-(rSi + lft- T2^)c 2 f N 2 f 

_(\ Lr\r w 3 _ iL — 4. i / S 

1 16 24 ^) ^f Iy f 192 A£ + 8 ^ 3 W ( 2 
-i- j- Fit Mr , nf M 2 R f M 3 



+ fifo + 4 1^/ + B 32 N f 



i 2 Nj + 4 3 Nj 



4 



Table 2: The terms in Eqs 

nf _ _39 
D 10 64 



.<20lforW,. 



, ■ - - --/c = 3. 
D ll 16 * T 16 5 



B 20 ~ 



* 32 6 T 64 ? 

7351 , 891 y ]e, [2177 ,351 / 1/:: 

2048 512 ? 3 J? 1.2048 512 ? 3 J? 

_ [221 I 81 r 1 z-3 , 531 M 81_ « 

[2048 512 ? 3 J? 2048? 1024? 

nf _ 101 £ _ [ _59_ j_ JL / 1 j?2 , .i f 3 il 

D 21 384 ? [ 1536 256 ? 3 J 5 "'"32? 512 

of _ _5_ c , _5_ t2 of [ 43 . 1 

-"22 ~~ ?R8 S 144 ? > D 33 ~~ 5184 """ 71 4 



Table 3: The /J^ M in Landau gauge (g = 0) for different N f and W c = 3. 



T6 1 

■2 

■e 



°22 
S 30 = 



!44 ? ' "33 _ L 5184 + 72 ^ 3 ] £ + 1728 

3791075 , 140271 y _ 246915 >- ] ^ 

98304 2048 ? 3 8192 ? 5 J ? 

r 459983 31887 y , 19035 y 1 

32768 4096 ? 3 4096 ? 5 I ? 
24561 _ 3681 
16384 



31887 y , 19035 y 1 
4096 ? 3 4096 ? 5 J? 



3681 r , 2475 y 

, 8192 ? 3 4096 ? 5 

, ' 11277 _ 6327 <■ , 7155 y 

8192 4096 ? 3 8192 ? 5 J? 
, 2727 , _513_ y _ _945_ r ] e5 _ : 
[4096 8192 ? 3 8192 "J? 1 

pf [ 203483 19821 y 2655 y ] y , [ 38057 
0/1^7/; 7n^fi S3 only S5 ? ' moo 



1569 
2048 
81 
256 



2619 rf> , 729 >7 
»i« - oivz— j- 16384 ? 16384 ? 

19821 y 2655 r\t A 38057 345 y 1 y-2 

2048 ? 3 2048 ? 5 J ? [ 12288 512 ? 3 J ? 

_ _39_ y 15_ «• 1 

128 ? 3 2048 ? 5 ? 



of _ [ 11 , 251 
^32 " ["36 + — 



, _21_ <■ _ M5_ y " y4 _ _5_13_ y5 , _63_ e6 

1024 ? 3 4096 ? 5 J ? 4096 ? 4096 ? 

251 / 1 t _ [i]9_ 3_ y ] *2 , _15_ pi 5_ y4 

384 ? 3 J? [3072 256 ? 3 J? 512? 128? 



where the bL all vanish in Landau gaugeQ The general expres- 
sions for the Z?f,„ are rather long and unilluminating, and hence 
not included here. The results for the phenomenologically most 
interesting case, N c - 3, however, are given in Tab. [2] 

The numerical values of the ySj VIM 's in Landau gauge are 
given, for N c - 3 and Nf = 0, 2 and 3, in Tab. [3] For the 
reader's convenience, we give also the numerical results for 
general £ and Nf = N c = 3, 



oMM 

Po 

nMM _ 

Pi - 

oMM _ 
Pi ~ 



oMM 
Pi 



(21a) 
(21b) 

(21c) 



(21d) 



= 2.25, 

4.0 - 0.421875^- 0.28125£ 2 + 0.140625£ 3 , 

20.9183 - 0.552182^ - 0.168435£ 2 

- 0.0187824£ 3 + 0.171387^ - 0.0791016£ 5 , 

160.771 + 10.5774<f- 2.46840£ 2 

- 0.145040£ 3 + 0.857841^ + 0.245698£ 5 

- 0.1 13708£ 6 + 0.0444946£ 7 , 



which results serve to illustrate the weakness of the £ depen- 
dence in the vicinity of Landau gauge. Note that, while the first 
coefficient, B™ M is gauge independent, and universal, the coeffi- 
cients beginning with B^ M are gauge dependent, as is typical of 
momentum subtraction schemes (as usual, the universal value 
of y6^ M is obtained only in Landau gauge). 



N f = 



N f 



N, 



oMM 

Po 


2.75 


2.4166 


2.25 




6.375 


4.79166 


4.00 




47.5075357 


29.1756592 


20.9183135 




392.7385102 


226.4690053 


160.7710385 



applications of QCD Green's functions within non-perturbative 
continuum approaches 111 2i 12211 based on DSEs or FRGEs, for 
example. The precise definition of the underlying renormalisa- 
tion scheme, the MiniMOM scheme, puts these approaches on a 
firmer ground, and should serve to resolve any previous misun- 
derstandings. We have already stressed its utility in providing 
a route to a lattice determination of a s requiring only a calcu- 
lation of two-point functions, which are relative easy to deter- 
mine with high precision in current simulations. Here we show, 
as an added bonus, that the MiniMOM coupling may provide 
a useful alternative to the MS coupling in certain phenomeno- 
logical applications. We do so by considering the expressions 
for the perturbative contribution to quantities relevant to a se- 
lection of phenomenological applications in the Nf = 3 regime, 
expanded in terms of either the MS or the MiniMOM coupling. 
With a^g 



af s (Q 2 )/n, a MM = af M (g 2 )/^, D u D 2 and D 3 



ZJ 3 + 5D,D 2 -5D 3 
(22) 



from Eqs. ( [15] ), and 

Ci = -D u C 2 = -D 2 + 2D\ and C 3 = - 
an observable, O, whose MS expansion is 

O = 1 + Ai atto + Ar a— + A 3 a— + 

1 MS MS J MS 

has an equivalent MiniMOM coupling expansion 

O = 1 h-AiAmm + [A 2 + CiAi]a 2 MM (23) 
+ [A 3 +2CiA 2 + C 2 Ai]a 3 MM 

+ [A 4 + 3dA 3 + (2C 2 + C\)A 2 + C 3 Ai] a 4 MM + . . . . 

In investigating the phenomenological utility of the Mini- 
MOM coupling, we will consider the N c = Nf = 3 case, for 
which Ci = -2.6875, C 2 = -2.38108495 and C 3 = 1.3815951. 
We show that, in some of the considered cases, use of the 
MiniMOM coupling significantly improves the apparent con- 
vergence of the relevant perturbative series, while, in other 
cases, it does not. Whether or not it is useful to employ the 
MiniMOM coupling is thus something to be decided on a case 
to case basis. 



4. Comparing perturbative expansions 

The definition of the running coupling in ([TJ has been widely 
used, and continues to be widely used, in phenomenological 



4 The N c = 3 Landau gauge version of these results were first presented in 
Refs. (^,01 and subsequently confirmed in Ref. 1191 . 



4.1. The Adler function of the vector/axial vector current cor- 
relators 

Our first example is the dimension D — contribution to 
the Adler function, D V/A . U (Q 2 ) = -Q 2 dn v/A]U (Q 2 )/dQ 2 , of 
the flavour ij vector (V) or axial vector (A) current scalar 
correlator, Tly/A-jj- At scales of phenomenological interest, 
Dv/A;ij(Q 2 )\ D _ is f ar ar, d away the dominant term on the OPE 



side of finite energy sum rules (FESRs) which have been stud- 
ied in the literature based on either electroproduction cross- 
sections or hadronic r decay data. The r-based FESRs are 
used in precision determinations of a ,(Mz ), the most recent of 
which are described in Refs. 1281 I29L 13011 . A combination of 
electroproduction- and r-based FESRs has also been used to in- 
vestigate the present discrepancy , 32 ] between the electro- 
production and t version of the / = 1 V spectral function 13311 . a 
discrepancy which prevents a clear decision as to whether or not 
the Standard Model (SM) prediction for (g - 2)^ is compatible 
with the current high-precision experimental result 13411 . 

Dv/A;ij\ D=0 is known to 5-loops l35ll and, for N c - Nf - 3, 
given in terms of a^, by 



An 2 D 



V/A;ij\ 



1+fl _ + 1 . 6398fl £_ 

+ 6.3710 a!- + 49.0757 at- + 

MS MS 



(24) 



The equivalent expansion in terms of omm is 



An l D 



V/A;ij\ 



1 +fl MM - 1-0477 at 



(25) 



• 4.8241 a 3 MM + 3.1257 a 4 MM + 



which displays significantly improved convergence at scales 
relevant to the phenomenological studies noted above (~ 2 — 
4GeV 2 ), even when one takes into account the increased size of 
a™ M {Q 2 ) as compared to a^ s (Q 2 ). This improved convergence 
will also be manifest in FESR studies which employ "contour 
improved perturbation theory" (CIPT) in their evaluations of 
the relevant weighted D — integrals |f| 



PS cases, apart from the overall (m, + rrij) 2 factors, are the dom- 
inant terms on the OPE side of S and PS FESRs and Borel sum 
rules (BSRs) which provide the most reliable sum rule determi- 
nations of m s and m u + mj Il37ll38ll39ll . Useful lower bounds on 
m s have also been obtained from the PS ij = us sum rules using 
a combination of the accur ately known K pole contribution and 
spectral positivity lf39l 4oL 4111 . A combination of FESRs and 
BSRs based on nps ;l(J , in addition, provides a determination, 
not just of m s , but also of the decay constants of the first two 
excited K resonances, and hence a useful, highly constrained 
model of the us PS spectral function, a model which, combined 
with the us S spectral function constructed in Refs. l38tl . al- 
lows the continuum 7 — spectral contributions to be sub- 
tracted from the experimental differential distribution in strange 
hadronic r decays. This turns out to be a crucial input to the 
hadronic r decay determination of \V US i il SI El III since 
the OPE representation of the 7 = contributions is extremely 
badly behaved at all kinematically accessible scales, prevent- 
ing one from employing FESRs based on the full experimental 
differential distribution [ 46] . 

The expansion of n^ ps . r in terms of a^ is known to five 
loops [39] and, for N c — Nf — 3, is given by 



3[( mi + mj)(Q 2 )] 2 



n 



8tt 2 Q 2 



11 

T 



1 + — fl MS + ( 26 ) 



+ 14. 17928 a* + 77.36834 at 



MS 



MS 



511.82848 a?— + 

MS 



where mi(Q 2 ) is the running quark mass in the MS scheme. Re- 
expressing the series in terms of %m yields 



4.2. The second derivative of the D—0 part of the scalar/ 
pseudoscalar correlator 

As our second example, we consider the subtraction- 
constant-independent second derivatives, n^p S . r (g 2 ), of the 
D — part of the scalar (S) and pseudoscalar (PS) correlators, 
Hs-,ij and Hps jj, formed from the divergences of the flavour ij 
V or A currents. These quantities, which are equal for the S and 



3[(m,+m J )(e 2 )] 2 



5 The CIPT prescription employs the truncated expansion of Eq. {25) point- 
by-point along the circle \Q 2 \ = so in the complex Q 2_ pl ane - An alternate 
approach to evaluating the weighted D = contour integrals is to use the 
"fixed order perturbation theory" (FOPT) prescription, in which the series is 
expanded, and truncated, using the running coupling at the same fixed scale, 
e.g., fi 2 = sq, for all points on the circle. In the FOPT scheme, large logs 
are unavoidable over some portion of the contour. Nonetheless, recent argu- 
ments 1 28], based on a model constrained by known features of the large order 
behaviour of the D = perturbative series for Dv/Ajj(Q 2 ), shows it is possi- 
ble that the truncated FOPT form might provide a more reliable representation 
of the resummed series than would the CIPT form. This is potentially rele- 
vant here because the improvement in the convergence of the integrated FOPT 
series achieved through the use of the MiniMOM couplant is, for commonly 
used weights, far less compelling than that achieved in the CIPT case. While a 
study of the 5-loop FOPT approximation to a range of weighted integrals of the 
model for the resummed series in Ref. 12811 shows that a good representation 
of the corresponding data integrals is not possible, in contrast to the situation 
when the 5-loop CIPT evaluation is employed [ 36], the analogous study has not 
yet been performed for the full resummed model, and, as a result, the prefer- 
ence for the CIPT approach (where the improvement due to the re-ordering of 
the series using a MM would be operative) is not yet conclusive. 



% PS ;ij(Q 2 )\ D=0 = 



8tt 2 Q 2 



1 + y«MM + (27) 



+ 4.325 12 a 



MM 



■ 7.57595 a 



MM 



71.99997 a 



MM 



which again displays significantly improved convergence. Such 
improved convergence is likely to allow a significant reduction 
in the errors on the determination of the light quark masses, and 
an improved version of the light quark mass bounds. 

4.3. The leading D — 2 contribution to the flavour-breaking 
ud-us, V+A, 7 = 0+1 correlator difference 

The flavour-breaking correlator difference, AI1 T = ^iy+A ud ~ 



n 



(0+1) 



where the superscript (0+1) denotes the sum of J 



and 1 contributions, is of interest in extracting |V HJ | (and/or 
m s ) from hadronic r decay data. The leading term in the OPE 
representation of AI1 T is the D - 2 mass-dependent perturbative 
contribution, AII T | , proportional to m 2 . FESRs based on the 
7 = 0+1 combination are employed because of the very bad 
behaviour of the OPE representation of the related integrated 
D - 2, 7 = contribution l46Tl . Even after the subtraction of 
7 = spectral contributions made possible by the us S and PS 
studies noted above 137113811 . the \V US \ (and/or m s ) extraction is 
complicated by the slow convergence of the D - 2, 7 = 0+1 

I OPE 

series at the correlator level. An T | D 2 is known to 4-loop order 



6 



and, for A^ c = Nf = 3, neglecting corrections of 0(m*/m%), has 
the form ||47ll 



2 .,ope _ 3m 2 (fi 2 ) 



2n 2 Q 2 



1 + 3^+ 19.93313 0^ 



+ 208.746 fl l_ + (2378) aL.+ 



(28) 



where the PMS/FAC estimate for the 5-loop coefficient, 
2378 [47], has been included for exploration purposes. Since 
fl Ms( m r) ~ 0.1, convergence is rather slow at the space-like 
point on the FESR contour. It would be very helpful in reduc- 
ing theoretical uncertainties in the determination of \V US \ by this 
method were the use of the alternate couplant, «mm, to signifi- 
cantly improve the convergence of the D - 2, 7 = 0+1 series. 
Recasting (l28l terms of omm, we find 



An T (g 2 )|°™ 



3<(g 2 ) 



2n 2 Q 2 



1 + -amm + 13.66230 a MM 



+ 96.04956 < M + (747.25429) a* + 



(29) 



Unfortunately, even at scales ~ m 2 = 3.16 GeV , where 
omm('«;) ~ 0.15, the decrease in the coefficient sizes for this 
alternate representation is roughly compensated for by the in- 
creased size of flMM as compared to a.y^. Thus, in this case, the 
MiniMOM coupling does not produce a useful re-ordering of 
the original series. 

5. Conclusions 

Our main intention here was to provide the precise definition 
together with a perturbative analysis of the MiniMOM scheme 
underlying the running coupling in ([]]) which has been widely 
used in non-perturbative studies of the infrared behaviour of 
QCD Green's functions and phenomenological applications 
since its introduction more than 12 years ago in Refs. [4]. 

The particular occasion for this probably overdue clarifica- 
tion is the relatively recent and promising effort to determine 
this coupling from lattice simulations. The high precision and 
reliable error estimates desirable for this project require the 
improved perturbative knowledge which our relation between 
the MiniMOM and the MS couplings at 4-loop order provides. 
When determining the MiniMOM coupling from the gluon and 
ghost propagators of lattice Landau gauge in Monte Carlo sim- 
ulations, with discretisation errors of 0(a 2 ), one has 



MM 



(<r) = 



g>) 

An 



Z L (q 2 ,a 2 )G 2 L (q 2 ,a 2 ) + 0(a 2 ), 



(30) 



as a — > 0, with g 2 (a) the bare coupling at the lattice cutoff scale 
1 /a, and Zi and Gl are the bare lattice gluon and ghost dressing 
functions, respectively (see ||5J for details). It turns out that 
a MU can k e determined quite accurately over a wide range of 
scales. Once the lattice scale is fixed, this translates into a high- 
precision determination of A MM . A MS can then be determined 
via 

A MM (DA 

= ex P — I ( 31 ) 



A MS 



2A> 



Table 4: A MM /A MS m Landau gauge for ^ c = 3 and different Nr. 



N f = 


Nf = 2 


N f = 3 


N f = 4 


Nf = 5 


1.8968 


1.8469 


1.8171 


1.7831 


1.7440 



where D\ is given by Eq. (1 1 5ab and, with our normalisation, 
for N c — N/ = 3 for example, /3o = 9/4. Numerical values 
for the N c = 3, Landau gauge version of this conversion factor 
are given in Tab. H A thorough up-to-date analysis using the 
available lattice data for Nf — 0,2 is given elsewhere @]. 

Finally, as we have demonstrated in the previous section, it 
turns out that the MiniMOM coupling is also well suited for 
use in perturbative analyses in the place of the commonly used 
MS coupling, since the relation between the two is known to 
4-loop order. We have compared the corresponding perturba- 
tive expansions of a selection of observables of phenomeno- 
logical interest and found that the MiniMOM coupling appears 
to provide an efficient re-ordering of the leading D — per- 
turbative contributions for the V/A and S/PS correlator cases. 
One should bear in mind, however, that this improvement is 
not universal, as the example of the D — 2 contribution to 
the flavour-breaking, 7 = + 1, V + A correlator difference 
shows, though this case does not create a dramatic problem ei- 
ther. Those examples where significant improvement is found, 
however, suggest it may be useful to consider q-^"^ as an alter- 
native expansion parameter for the perturbative series relevant 
to other observables, such as those relevant to the decays of 
heavy quarkonia and heavy quark physics, as well. 
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